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Abstract 


The envelope theory, also known as the auxiliary field method, is a simple technique to compute 
approximate solutions of Hamiltonians for N identical particles in D-dimension. The accuracy of 
this method is tested by computing the ground state of N identical bosons for various systems. A 
method is proposed to improve the quality of the approximations by modifying the characteristic 
global quantum number of the method. 
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I. INTRODUCTION 


The envelope theory (ET), aslo known as the auxiliary field method, is a tec 
obtain approximate solutions of A^-body Hamiltonians in quantum mechanics [l|, 


mique to 
It has 


been recently extended to treat systems with arbitrary kinematics in D-dimension [^. The 
basic idea is to replace the Hamiltonian H under study by an auxiliary Hamiltonian H which 
is solvable, the eigenvalues of H being optimized to be as close as possible to those of H. 
The method is easy to implement since it reduces to find the solution of a transcendental 
equation. In the most favourable cases, the approximate eigenvalue is an analytical lower or 
upper bound. In less favourable situations, a nonvariational numerical approximation can 
be computed. This is often interesting for A^-body problems which are always difficult and 
very heavy to solve numerically. 

In a recent paper accurate ground state energies and correlation functions of bosons 
interacting via various potentials are calculated using Explicitly Correlated Gaussian (ECG) 
basis [^. Up to 8 particles are considered, allowing a plenty of numerical results to test the 
accuracy of the ET for eigenvalues and eigenvectors (the quality of eigenvectors for the ET 
has only been tested for A^ = 2 Q|). It is worth mentioning that the purpose of the ET 
is not to compete with methods such as the ECG. The main goal of the ET is to produce 
without great pain reliable estimations for the energy of a A^-body system. If a lower or an 
upper bound can be produced, this must be sufficient for some applications or can be used 
as a test for heavy numerical computations. In the peculiar situations where an analytical 
expression is obtained, the result can give valuable insights about the system. 

In several cases, it has been shown that the approximation given by the ET can be 
largely improved by modifying the structure of the characteristic global quantum number 


B 


of the method 17|, 18|. The problem is that the kind of modification must be guessed for 
each system, and can be found precisely only by comparing with several exact (or accurate 
estimation of) eigenvalues. Clearly, a better defined procedure would be preferable. A 
universal effective quantum number for centrally symmetric 2-body systems is proposed in 
Qj. We will show that this notion can be used to improve the ET. Let us note that the 
variational character of the approximation cannot be guaranteed when the global quantum 
number is modified. 

The paper is organized as follows. The ET method is recalled in Sec. [TTl where a modifica- 
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tion of the global quantum number is proposed to improve the quality of the approximation. 
The accuracy of the method is tested with various systems in Sec. IIIIl Some concluding 
remarks are hnally given. 


II. THE ENVELOPE THEORY 


Here we consider the V-body Hamiltonian, in a D dimensional space {D > 2), for 
identical particles with a kinetic energy T, interacting via the one-body U and two-body V 
interactions (h = c = 1) 


N N N 

H = J2nM) + j2u{\ri-R\)+ 

i<j=l 


( 1 ) 


2=1 2=1 


where P* = 0 ^ind R = ^ Sj=i centre of mass position. It is shown in 

that an approximate eigenvalue E is given by the following set of equations for a completely 
(anti)symmetrized state: 

To 


roPo = Q, 

Np„r(p„) = r„f/' (^) + ^r„V’ (^) , 


\VCn J 

where A'{x) = dA/dx, Cm = N{N — l)/2 is the number of particle pairs, and 

N-l ^ 


Q — -|- /j) -t- {N — 1)-— 


( 2 ) 

(3) 

(4) 

(5) 


2=1 


is a global quantum number (corresponding to V — 1 identical harmonic oscillators). Fol¬ 
lowing the forms of T, U and V, the approximate value E can have a variational character. 

n 

The corresponding approximate eigenstate is given by [8| 

N-l 

V’ = 0cm(-R) Yl ^ndiiKXi), (6) 


2=1 


where cCj) is a H-dimensional harmonic oscillator wavefunction [^, depending on the 

Jacobi coordinate Xi j^, and decreasing asymptotically like (^^he magnetic quantum 

number is omitted). The function 4>cm{R), which can be an oscillator state, describes the 
centre of mass motion {xm = R)- The scale parameters A* are given by 


A* = 


i N 


—NQ— = 

i + 1 ro \l i + 1 Q 


Po- 


(7) 
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A state (j6]) has neither a dehned total angular momentum nor a good symmetry, but its is 
characterized by a parity (—By combining such states with the same value of 
Q, it is generally possible to build a physical state with good quantum numbers and good 
symmetry properties, but the task can be technically very complicated, even for N = 3 11 |. 


It is possible to compute mean values for the power of interparticle distance (r^) with 
the pair correlation function dehned by 

N 


C{r) = 


'ip 


i<j=l 


r) 




( 8 ) 


N{N -1) 

The value at origin, 6, of the radial part of this density is also interesting to compute. For 
the bosonic ground state (n, = k = 0, Vi), which is completely symmetrical, we have 

2Af 


C{r) = 


and 5 = 


( 9 ) 


In 


“ T{D/2)' 

9(], it is shown that an effective quantum number q that determines with high accuracy 
the level ordering of centrally symmetric 2-body systems has the following structure)? ] 

2 


, = ,^(n + i)+! + k_ 


( 10 ) 


where the number 0 depends on the system. For instance, exact results are obtained for a 
harmonic oscillator with </> = 2, and for a Coulomb system with 0=1. In the dominantly 
orbital state (DOS) method [12|, ll3|, the approximate solutions are found by quantizing the 
radial motion around a semiclassical solution for a circular motion. The same separation in 
a radial quantum number n + ^ and an orbital quantum number I + naturally occurs 
in this method. Its extension to 3-body systems seems to preserve the same separation 13 |. 


Moreover, it has been shown, in the framework of the FT, that relations exist linking the 
energy of a A^-body system to the energy of the corresponding 2-body system by a rescaling 
of the global quantum numbers IJ]- So, it is tempting to try to modify the global quantum 
number Q as 

D + (b-2 


N-l 

+ h) + (N — 1 )- 


( 11 ) 


2=1 


which partly breaks the strong degeneracy of Q. The parameter 0 can be determined by 
theoretical considerations, or by a £t on a single known accurate solution. In the following, 
Q will be systematically replaced by Q,/,. The genuine FT, with its possible variational 
solutions, is then recovered with 0 = 2. For other values of 0, the variational character of 
the solution cannot be guaranteed. 
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III. RESULTS 


In this section, we consider various systems of N identical particles with the same mass 
m. The ET formalism developed is valid for bosons or fermions with any values of D, but 
we will show only numerical results for bosons in the D = 3 space. For the three hrst 
nonrelativistic systems taken from and named as in only the ground state is computed. 
The fourth case is an ultrarelativistic 3-body system studied in [TJ , for which excited states 
are also available. 


A. Weakly interacting bosons 


In the hrst case, the nonrelativistic particles interact via a soft Gaussian potential 

T(P) = 1^, U(s) = 0, V(r) = 

Im 

The resolution of the system (EH) gives the following upper bound for the energy E 

ro = N^I\N R, 

E = 

2 

1 


N{N-1) + 


Wo{Yy 


withr = 


iVi/2 (iV - l)Ry/27^' 


( 12 ) 

(13) 

(14) 

(15) 


The multivalued Lambert function W(z) is the inverse function of z 15|. ITo(^) is the 
branch dehned for z > — 1/e. 

In order to test the quality of the approximation with accurate data, we compare with 


the results given in 


for the parameters 1/m = 43.281307 


a.u., 


K, Vo = 1.227 K and 


R = 10.03 a.u. Results for the 


round state energies E are presented in Fig. [T] (the exact 


result for TV = 20 is taken from lOj). The agreement is quite good for the genuine ET, but 


irrelevant values for the energies (complex or positive real numbers) are computed when the 
binding is very weak. This phenomenon has already been observed in for N = 2. Note 


that excited values of the energies are available up to A = 6 in 


18| . but the ET method can 


only produce one relevant value of the energy for N = 6. So no tendency can be shown as a 
function of N. A better agreement with the exact values for A > 5 is obtained for 0 = 1.82, 
but the improvement is not spectacular. 

Mean values (r) are presented with 6 in Fig. |2l The global agreement is good, but the 
value 0 = 1.82 does not necessarily produce an improvement. Note that the values of 6 are 
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Figure 1. Binding energies, —E, of weakly-interacting bosons: exact results (circle) 4, Q, ET 
results for <p = 2 (diamond), ET results for cf) = 1.82 (dashed line). Left: log-log plot, right: linear 
plot. 


not so well reproduced. This is typical for a given value of a density. Better results are 
obtained for mean values, for which an integration is performed with the density over the 
whole domain of distances. 
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Figure 2. (r) (left) and 5 (right) for weakly-interacting bosons: exact results (circle) ET results 
for (j) = 2 (diamond), ET results for cj) = 1.82 (dashed line). 


B. Self-gravitating bosons 


In this case, we consider nonrelativistic particle interacting via a Coulomb-like potential 


T{p) 


P 

2m 


, U{s) 


0, V{r) 


9 


r 


(16) 


6 











For a self-gravitating system, g = m^G where G is the gravitational constant. The ET 
method gives the following upper bound 


23/2 g2 
^0 iVV2(iV- 1)3/2^^’ 


For iV = 2, the exact result is recovered with <p 
the ground state bosonic energy is 


1 10|. Moreover, for D 


(17) 

(18) 

3 and 0 = 1, 


Egs ^ -0.0625 N^{N - 1) m /. 
This is comparable with the lower bound obtained in 


(19) 


Egs >-0.0593 l)m/. (20) 

As in the previous case, the comparison is performed with j4|, for the parameters m = 1 
and g = 1. Results for the ground state energies E are presented in Fig. |3]for 0 = 2 and 
0 = 1. If the global tendency is good for the genuine ET, as indicated by the log-log plot, 
a spectacular improvement is reached with 0 = 1 (a better global agreement with the exact 
values for 2 < iV < 8 is obtained for 0 = 1.11). 



Figure 3. Binding energies, —E, of self-gravitating bosons: exact results (circle) ET results for 
0 = 2 (diamond), ET results for 0 = 1 (dashed line). Left: log-log plot, right: linear plot. 


Mean values (r) are presented with 6 in Fig. m Again, for the genuine ET, values of (r) 
are globally better reproduced than values of 6. The strong improvement obtained for the 
energy is not repeated for these observables. 
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Figure 4. (r) (left) and S (right) for self-gravitating bosons: exact results (circle) 
for (j) = 2 (diamond), ET results for (f> = 1 (dashed line). 



ET results 


C. Confined bosons 


Now, we consider particles confined by a harmonic oscillator potential with a pairwise 
repnlsive Conlomb interaction 

.2 1 


T{P) = U{s) = -mu^s\ V{r) = L 


The ET method gives the following lower bonnd 

7V5/6(iV- 1)V2 . g y/3 


Tq = 

E = 


25/6 

7V2/3(/V- 1) 


mu 
2 „2^V3 


) G_(F), 


withy = 


25/3 

216/3 


(m g^) (G_(y)2 + G_(y)-i) 


u 


2/3 


Ql 


( 21 ) 

( 22 ) 

(23) 

(24) 


3 iVy5(_/V —1)2 

G±{Y) is the only positive root of the qnartic eqnation 4a;^ ± Sx — 3y = 0 with y > 0 
[? ]. One can check that E = uQ when S' = 0, which is the exact solntion in this case. 

Again, the comparison is performed with j4|, for the parameters m = 1, cu = 0.5 and g = 1- 
The Hamiltonian stndied in this reference is dehned with Xlili '^1, — + NR, 

that is to say with the centre of mass motion rnled by a harmonic oscillator, with a mass 
Nm and a freqnency u. So, to compare the resnlts of this paper with onrs, we mnst add the 
valne |a; (centre of mass gronnd state energy) to the internal energy of the gronnd state of 
onr Hamiltonian (ET]). Resnlts for the gronnd state energies E are presented in Fig. [5l The 
agreement is qnite good for the gennine ET. With the valnes chosen for the parameters, an 


analytical solntion can be fonnd for N = 2 


20| . This exact resnlt can be reprodnced by onr 
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formula with 0 = 2.58. A better global agreement is obtained for this value of 0, but the 
improvement is very weak. 



Figure 5. Energies, E, of confined bosons: exact results (circle) ET results for </> = 2 (diamond), 
ET results for (j) = 2.58 (dashed line). Left: log-log plot, right: linear plot. 


Mean values (r) are presented with 5 in Fig. |6l The agreement of the genuine ET results 
is just reasonable, with again the values of (r) better reproduced than the values of 5. In 
this case, the quality of the observables is deteriorated for 0 = 2.58. 



Eigure 6. (r) (left) and 5 (right) for confined bosons: 
(diamond), ET results for = 2.58 (dashed line). 


exact results (circle) , ET results for (^ = 2 
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D. Large-A^ baryons 


Light baryons in a particular large-iV limit {N is the number of colours) can be studied 


with the following Hamiltonian 




T(p) =p, U{s) = As, V{r) = - 


9 


(25) 


The baryons are then composed of N quarks u or d, with m = 0, moving ultrarelativistically. 
Quarks are fermions, but the colour part of the wave-function is completely antisymmetrical. 
So, they behave like bosons for the rest of the wave-function. With the ET, the following 
upper bound is obtained for the mass E of the baryons 


^0 = ^ \Inq^ 




3/2 


9 , 


E = VIX\ Ni 


(«^) 


3/2 


9- 


(26) 


(27) 


In order to test the relevance of for the ground state and some excited states, we 
consider the 3-body system studied in [t]. For the genuine ET (0 = 2), the mean relative 
error A computed on the 16 eigenmasses listed in Table [I] is 15.1%. With the value 0 = y/2 
predicted by the DOS method ioi g = 0 Q, A reduces to 4.4%. With 0 = 1.35, the ET 
yields the exact value for the ground states and A drops to 3.1%. Finally, with 0 = 1.23, 
the minimal value of 2.4% is reached for A. 


IV. CONCLUSION 


The envelope theory is a powerful method to compute eigenvalues for quite general N- 
body systems with identical particles in D dimensions . The method is easy to implement 
since it reduces to hnd the solution of a transcendental equation. Its interest is not to 
produce accurate results, but to yield rapidly reliable estimation of energies and observables. 
Here, the method is tested with the ground states of three nonrelativistic systems, weakD 
interacting bosons, self-gravitating bosons, and conhned bosons, accurately computed in |4| 
up to A = 8. A fourth case is an ultrarelativistic system of quarks (boson-like) studied 
in [^. For all these cases, an analytical variational bound is computed. Considering the 
simplicity of the method, quite good results can be obtained for energies and observables in 
all the systems considered. 
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A universal effective quantum number for centrally symmetric 2-body systems is proposed 
in jo I. As it seems relevant to extend this notion to 3-body systems j^|, the characteristic 
global quantum number of the envelope theory is modified accordingly. Only one parameter 
must be hxed either by theoretical considerations or by a fit on a single known accurate 
solution. The variational character of the energies cannot then be guaranteed, but the 
improvement is generally good, even for the excited states. Unfortunately, the behaviour of 
observables is not really predictable. 


Since A^-body problems are always difficult and heavy to solve accurately, the envelope 
theory can be used as a guide for the study of these complicated systems. It could be 
interesting to test the method with other systems in various dimensions, for different values 
of N, and not only for the ground state. In particular, it is not certain that the modification 
of the global quantum number proposed in this work can always lead to improvement of the 
energies. 
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basis 


Table I. Some eieenmasses in GeV of the Hamiltonian ( 1251) for D = 3, N = 3, X = 0.2 GeV and 

n 

g = with as = 0.4 [T]. Accurate results obtained from an expansion in a harmonic oscillator 


22l ] are given in the column “Exact”. Results computed with the ET for four values of cj) are 


listed in the four last columns, with the associated mean relative errors A. The sums ni + n 2 and 
li + I 2 are the quantum numbers of the main component of the corresponding genuine eigenstate 
in the harmonic oscillator basis. These numbers are used for the computation of 


ET 


ni + 

n2 h 

+ I 2 Exact 

(/> = 2 

(l) = y^ (j) = 1.35 (j) = 1.23 

0 

0 

2.128 

2.468 

2.165 

2.128 

2.060 

0 

1 

2.606 

2.914 

2.662 

2.633 

2.578 

1 

0 

2.739 

3.300 

2.842 

2.788 

2.682 

0 

2 

2.959 

3.300 

3.080 

3.055 

3.007 

1 

1 

3.125 

3.646 

3.237 

3.189 

3.098 

0 

3 

3.299 

3.646 

3.448 

3.425 

3.383 

2 

0 

3.260 

3.961 

3.387 

3.318 

3.186 

1 

2 

3.422 

3.961 

3.589 

3.546 

3.463 

0 

4 

3.581 

3.961 

3.780 

3.759 

3.721 

2 

1 

3.584 

4.253 

3.725 

3.662 

3.542 

1 

3 

3.716 

4.253 

3.909 

3.869 

3.794 

0 

5 

3.861 

4.253 

4.085 

4.066 

4.030 

3 

0 

3.721 

4.527 

3.856 

3.775 

3.619 

2 

2 

3.838 

4.527 

4.034 

3.976 

3.866 

1 

4 

3.966 

4.527 

4.205 

4.168 

4.098 

0 

6 

4.103 

4.527 

4.369 

4.351 

4.318 

A 



15.1% 

4.4% 

3.1% 

2.4% 
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